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Abstract. This paper develops a geometric model for coupled two-state quan- 
tum systems (qubits) using geometric (aka Clifford) algebra. It begins by 
showing how Euclidean spinors can be interpreted as entities in the geometric 
algebra of a Euclidean vector space. This algebra is then lifted to Minkowski 
space-time and its associated geometric algebra, and the insights this provides 
into how density operators and entanglement behave under Lorentz transfor- 
mations are discussed. The direct sum of multiple copies of space-time induces 
a tensor product structure on the associated algebra, in which a suitable quo- 
tient is isomorphic to the matrix algebra conventionally used in multi-qubit 
quantum mechanics. Finally, the utility of geometric algebra in understanding 
both unitary and nonunitary quantum operations is demonstrated on several 
examples of interest in quantum information processing. 



1. Introduction 

Quantum mechanics attaches physical significance to representations of the ro- 
tation group which differ substantially from those studied in classical geometry. 
Much of the mystery surrounding it is due to this fact. The enormous interest re- 



cently generated by proposals to build a quantum computer |Llo95, EJ96, 3te98, 
WC98, Bro9g| , BDOCl| has focussed attention on the simplest possible quantum 



system: a two-state system or qubit. Our understanding of qubits is based on two 
distinct geometric models of their states and transformations: 

• A complex projective line under the action of SU(2). 

• A Euclidean unit 2-sphere under the action of SO (3). 

The first is used almost exclusively in fundamental quantum physics, while the 
second ("classical") model is used extensively in certain applications, e.g. nuclear 



magnetic resonance (NMR) spectroscopy |HSTCOO|. In particular, in quantum 
computing a qubit represents a binary or 1 as its state corresponds to one of a pair 
of conjugate (orthogonal) projective points | 0) or | 1). These in turn correspond 
to a pair of diametrical points on the unit sphere, which determine the alignment 
of the qubit with or against the corresponding axis of quantization. 
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Formally, these two models are related by stereographic projection of the Rie- 
mann (unit) sphere onto the Argand plane, the points of which are the ratios of the 
homogeneous coordinates of points on the projective line (see e.g. [ Alt86, FH81|). 
While this elegant construction describes the mapping between the two representa- 
tions in a geometric fashion, it does not unite them in a single mathematical struc- 
ture. This paper provides an informal account of how this is done by geometric (aka 
Clifford) algebra; in addition, it describes an extension of this formalism to multi- 
qubit systems, and shows that it provides a concise and lucid means of describing the 
operations of quantum information processing | SCH98| , HSTCOC]. Significantly, 
this extension is most naturally derived via the geometric algebra of Minkowski 



space-time [DLG93|, which has also been shown to be an cfhcicnt formalism within 
which to study a very wide range of problems in classical |Hes99 , Jan89|, rclativis- 



tic 



Hes66| Bay96] and fundamental quantum DLG+96| physics. More complet e 



and rigorous accounts may be found in these references, and in [ HCSTOO , HavOlj 



2. Euclidean Geometry and Spinors 

Let 9^3 be a three-dimensional Euclidean vector space whose inner product is 
denoted hy {a,b) i-^ a ■ b. The Clifford or geometric algebra of 9^3 is the associative 
algebra generated by Dl^ over M such that a^ = ||a|p = a- a for all a S 9^3. This 
algebra will be referred to in the following as the Pauli algebra, and denoted by ^^3. 
The interesting thing about this algebra is its geometric interpretation, which will 
now be described. 

To begin, note that every nonzero vector a G 9I3 has an inverse a/||a|p. In 
addition, a simple application of the law of cosines shows that the inner product 
of a with any other vector b G SH3 is given by the symmetric part of the geometric 
product: 



(2.1) 



i(ab -|- ba) 



bf 



The antisymmetric part, by way of contrast, is called the outer product, and denoted 
by (a, b) i~> a A b = {ab — ba)/2. Since the outer product of two vectors a l\b 
is invariant under inversion in the origin, it cannot itself be a vector. The space 
(a A b I a, b G IHs) therefore carries an inequivalent representation of the orthogonal 
group 0(3), and its elements are accordingly called bivectors. These are most 
naturally interpreted as oriented plane segments, instead of oriented line segments 
like vectors in £H3 . If we similarly define the outer product of a vector with a 
bivector and require it to be associative, i.e. 



(2.2) 



a A {b Ac) = i (abc - cba) = {a Ab) Ac 



(a, b,c € 9I3), then it can be shown via straightforward though somewhat lengthy 
calculations that this product of three vectors is totally antisymmetric, meaning 
that the outer product generates the well-known exterior algebra /\ 9I3 (cf. |Hes99, 
Rie58|). The outer product of three vectors is called a trivector, and (since it 



changes sign under inversion) is most appropriately interpreted as an oriented space 
segment or volume element. 

The general properties of inner and outer products in the geometric algebras of 
arbitrary metric vector spaces can be worked out along these lines in a coordinate- 
free fashion [HS84|. The remainder of this section will focus on how the Pauli 
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algebra is used to describe the quantum mechanics of qubits. In this apphcation it 
is more common to work with a fixed orthonormal basis (jyi^^ (jy ^ ct-^ e 9I3. Quantum 
mechanics, however, views these basis vectors in a very different way from that 
taken above, in that they are regarded as operators on a two-dimensional Hilbert 
space i3 w (see e.g. [Sak94|). These operators, in turn, are usually identified 
with the Pauli matrices 



(2.3) 



I 

1 



-I 

1 





-I 



where i is an imaginary unit (r^ — —I), the underline signifies that the associated 
symbol is a matrix, and throughout this paper the symbol "<-!■" should be read as "is 
represented by" or "is equivalent to" . The connection between the two viewpoints 
lies the fact that these matrices satisfy the defining relations of the abstract Pauli 
algebra Q3, namely 



(2.4) 



1 



and hence constitute a faithful matrix representation of it. This shows, in particular, 
that is 8-dimensional as a real linear space.^ 

In most physical situations, these operators (times h) represent measurements 
of the intrinsic angular momentum of the qubits, and hence are regarded as genera- 
tors of rotations in the Lie algebra so(3) over C satisfying the commutator relation 



(2.5) 



and its cyclic permutations. In terms of geometric algebra, the left-hand side is 
just the outer product of the vectors. The right-hand side is somewhat harder to 
interpret, because the Pauli algebra is defined over the real numbers. The trick 
is to observe that, in terms of the matrix representation, g_^g_yg_^ = «1. Thus by 



interpreting the abstract imaginary i as the trivector i 
the angular momentum relations become a triviality: 

n2 



(Tx (Ty CTz — C7"x A CTy A (T^ 



(2.6) 



More generally, the vector cross product is related to the outer product by 



(2.7) 



axb = —^{ab—ba) = — <,(a A b) 



from which it may be seen that multiplication by the unit trivector l maps vectors 
to orthogonal bivectors and vice versa. Since they span a one-dimensional space but 
change sign under inversion in the origin, trivectors can also be regarded as pseudo- 
scalars. Perhaps the most important thing which geometric algebra contributes 
to physics are geometric interpretations for the imaginary units which it otherwise 
uses blindly. 

If we denote the induced bivector basis by 



(2.8) 



(Ty A 



(T, A (Tx 



K 



<Tx A cr, 



In the quantum mechanics Uterature, the notation a - is often used for f^M^/i ■ Because 
a is geometrically just a vector in SH3 {not a matrix for it in a basis-dependent representation 
of 63 ), this is an abuse of the dot- notation for the Euclidean inner product, which is otherwise 
perhaps the most consistently used notation in all of science. This abuse of notation will not be 
perpetrated in this paper. 
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it is readily seen that these basis bivectors Ukewise square to —1. On muhiplying 
the angular momentum generating relations through by — 1 = t^, we obtain 

(2.9) JI ^ K , IK ^ J , KJ = I , and KJI = 1 . 

This shows that these basis bivectors generate a subalgebra of Gs isomorphic to 
Hamilton's quaternions | Alt86 , Alt89|, which is also known as the even subalgebra 

(since it is generated by the products of even numbers of vectors) . It is well- 
known that the quaternions' multiplicative group is R* © SU(2), which implies that 
the even subalgebra should be closely related to rotations. This relationship will 
now be worked out explicitly. 

Consider the result of conjugating a vector a; by a vector a, i.e. 

(2.10) axa^^ = a{x,, + x^)a^^ = aa^^x,, — aa^^x^^ , 

where we have split x — {x ■ a^^ + x /\ a^^)a = a;,, + a;^ into its parts parallel and 
perpendicular to a. This shows that axa^^ is the reflection of x in the plane 
orthogonal to a. From the well-known fact that the composition of two reflections 
is a rotation by twice the lessor angle between their planes and about these planes' 
line of intersection, it follows that conjugating a vector by an element of the even 
subalgebra just rotates it accordingly: 

, , , ,-1 -1-1 baxab 

(2.11) (ba)x(ba) = baxa b = j. — ,,„,,, ,,„ 

Let u = a/||a||, v = and R = vu be the corresponding unit quaternion. 

Then R = cos(6'/2) — trsin(6'/2) where cos(0/2) = u ■ v and tr = u /\ v/\\u A v\\. 
Moreover, the inverse (vu) ^ is now simply the reverse uv = {vu)\ which in 
turn corresponds to the conjugate quaternion R^ = cos(0/2) -I- Lrsm{9/2). This 
reversal operation on G-^ extends to a well-defined anti-automorphism of Gs , which 
corresponds to Hermitian conjugation in its representation by Pauli matrices. On 
splitting X into its parts parallel x„ and perpendicular Xj^ to r as above, the rotation 
may now be written as RxR} — R{x„ + x^)R^ = x„ + x^{R^)'^ 

= a;,, + a;^(cos2(6i/2) - sin2(6l/2) + 2<.r cos(6l/2) sin(6i/2)) 

= a;,! + a;^(cos(t^) + tr sin(6')) , 

and so may be viewed as multiplication of x^ by the "complex number" cos(6') -|- 
Lr s'm{0) in the Argand plane defined by the bivector tr. 

By collecting even and odd powers in its Taylor series, it may be seen that any 
unit quaternion can be written as the exponential of a bivector orthogonal to the 
axis of rotation r: 

(2.13) e-'"-"/^ = i_,rf - i(f )' + ••• = cos((?/2) - tr sin(0/2) 

This is formally analogous to a complex exponential, and is also in accord with 
our previous observation that the space of bivectors is isomorphic to the Lie alge- 
bra su(2) w so(3) under the commutator product. The pair [cos(6'/2); sin(6'/2)r] 
are known as Euler-Rodrigues parameters for the rotation group S0(3); since 
[— cos(0/2); — sm(0 /2)r\ determines the same rotation, this parametrization is two- 
to-one. A one-to-one parametrization is obtained from the outer exponential^ i.e. 

(2.14) A"''''^ = 1 - trr - ir A rr^ H = 1 - trr (since r A r = 0) . 

The squared norm of this outer exponential is 1 -I- , so that the normalized outer 
exponential equals the usual exponential if we set r = tan(6'/2). Because t = 
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tan(6'/2)r is the four-dimensional stereographic projection of [cos(6'); sin(0)r] from 
6* = TT, it has been called the stereographic parameter for S0(3). Note, however, 
that this parametrization does not include rotations by tt. 

Another two-to-one parametrization of rotations is given by the Cayley- Klein 
parameters [ipi\'ip2] G C^, where 



(2.15) 



V-i = cos(6'/2)-^sin(6'/2)(r•c^,) , 

'4>2 = sin(e'/2) (r • (Tx) +isin(0/2) (r • cTy) 



The corresponding SU(2) matrix is simply 



(2.16) 



■^1 

V'2 



-r2 



It follows that the complex column vector | V") = [V'i;'02] itself transforms under 
left-multiplication with matrices in SU(2), which is commonly described in quantum 
mechanics by calling it a spinor. In particular, the spinors | ) = [1; 0] and | 1 ) = 
[0; 1] are those commonly used in quantum computing to store binary information. 
Since the Cayley-Klein parameters uniquely determine the SU(2) matrix, however, 
we can just as well regard spinors as entities in SU(2), e.g. | ) <-!■ 1 and | 1 ) <-!■ 
The usual action of SU(2) on spinors then becomes the left-regular action of SU(2) 
on itself. 

The representation of SU(2) used above depends upon the choice of coordinate 
system: Changing to a different the coordinate system gives a different (though 
equivalent) representation. Recalling that SU(2) is isomorphic to the multiplicative 
group of unit elements (quaternions) in the even subalgebra , a coordinate-free or 
geometric interpretation of spinors is obtained by regarding them as elements of 
itself. This interpretation of spinors as entities in ordinary Euclidean geometry was 
first pointed out by Hestenes over thirty years ago [ Hes66 |, but physicists persist 
in putting operators and operands into separate spaces, and in working with a 
matrix representation instead of directly with the geometric entities themselves. 
The perceived nonintuitive nature of quantum mechanics is due in large part to the 
resulting confusion over the geometric meaning of the objects with which it deals, 
which is spelled out explicitly in geometric algebra. 

As another example, consider how the density operator of an "ensemble" of 
qubits can be interpreted in geometric algebra. This operator p is usually defined 
via a matrix representation as p = | -0 ) ( -0 | , where the over line denotes the average 
over the ensemble. As first observed by von Neumann, this matrix contains all 
the information needed to compute the ensemble average expectation values of the 
qubit observables, since 



(2.17) 



= tr(^^|7A)(0|) - tr(^^|V')(V'l) = tr(^^p) 



(/i G {x, y,z}). To translate this into geometric algebra, we set the second column 
o{ ^ to zero by right-multiplying it by the idempotent matrix = {X + g;^)/2, 
i.e. 



(2.18) 





-r2 




'1 0' 




0" 













i,2 0_ 
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This corresponds to projecting ^ G onto a left-ideal in Q3 , and allows the 
dyadic product in Eq. ( 2.17 ) to be written as: 



(2.19) 



>1 


6 








_i>2 














= (*^+)(*^- 



Thus the interpretation of the density operator in geometric algebra is 

(2.20) p = (*£;+)(*£;+)^ = " K"^ 



(cf. |SLD99|). The vector part p = (p)-^ = 'S'tTz'S't is called the polarization 
vector (in optics, its components are known as the Stokes parameters [BBDH93|, 



while in NMR it is known as the Block vector after the pioneer of NMR who 
rediscovered it |Blo46|). Its length is < 1 with equality if and only if all 
members of the ensemble are in the same state In this case the ensemble is said 
to be in a pure state, and the density operator is itself an idempotent (1 +p)/2, 
where p = '^cr^'^K For an ensemble in a general mixed state, the length of the 
ensemble-average polarization vector measures the degree of alignment among the 
(unit length) polarization vectors of the individual members of the ensemble, and 
is called the polarization of the ensemble. 

In many physical situations there is a natural reference direction; for example, 
in NMR computing the qubits are spin 1/2 atomic nuclei whose intrinsic magnetic 



dipoles have been polarized by the application of a strong magnetic field [ HCSTOO . 
From a geometric perspective, however, the density operator is just the sum of a 
scalar and a vector, which for a pure state is related to the corresponding "spinor" 
by rotation of a fixed reference vector (conventionally taken to be cr^ as above) 
by 'S'. Since the trace in the standard matrix representation is simply twice the 
scalar part (_)q of the corresponding expression in geometric algebra, the ensemble- 
average expectation value 

= cr • (*crz*^) = 



(2.21) 



is just the component of the polarization vector along the /i-th axis. Unlike the 
strong measurements usually considered in quantum texts, where measurement of 
07, yields one of the random outcomes ±1 with probabilities {Izter^ ■p)/2 and leaves 
the system in the corresponding state p = ±cr^ , weak measurements of ensemble- 
average expectation values can be made with only negligible perturbations to the 
ensemble as a whole [ Per93 1 . This is in fact how quantum mechanical systems are 
usually manifest at the macroscopic level! 

To see how all this relates to conventional wisdom, observe that the polarization 
vector of a pure state may be written in terms of the Cayley-Klein parameters as 

(2.22) p - 23fi(^i>2)crx + 2^irM<Ty + (|Vi|' - |V'2|')crz . 

Its stereographic projection from — ctz onto the cr^cry plane is therefore 

23fi(7/;i*^2)Tx + 2^{rM(Ty 



(2.23) 



l+|^lP-|^2P 



Multiplying by cr^ and simplifying the denominator using |^i 



1^2 



1 yields 



(2.24) 
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where K = cr^cry is a square- root of —1. This is the same as the ratio ?A2/V'i save 
for the use of i instead of K as the imaginary unit, which explains formally why 
S0(3) acts on the polarization vector in the same way that SU(2) acts on the ratio 



of the Cayley-Klein parameters |Alt86, FH81 



3. Space-Time Geometry and Multiparticle Spinors 



The above interpretations apply only to single qubits (or to ensembles consist- 
ing of noninteracting and identical qubits). Extending them to systems of interact- 
ing and distinguishable qubits may be done in a physically significant fashion by 
considering the geometric algebra of space-time (or Minkowski space) ^1^3 ■ This 
algebra, known as the Dirac algebra and denoted by may be defined by the 
generating relations among an orthonormal basis analogous to Eq. (2.4): 



(3.1) 



7? 



1 



7, 



-1 in e {x,y,z}) , 



7^7,. -7i.7^ (m,*^ e {t,x,y,z},/x 7^ i/) 



The corresponding geometric algebra separates into five inequivalent representa- 
tions under the action of the full Lorentz group 0(1, 3), i.e. 



(1) 



(3.2) 



(7^) 

(7t7x7y7z> 



(scalars, 1-dimensional) 
(vectors, 4-dimensional) 
(bivectors, 6-dimensional) 
(trivectors, 4-dimensional) 
(pseudo-scalars, 1-dimensional) , 



where ^,v^ri ^ {t, x, y, z} with \i ^ v ^ j] ^ [i^ for a total dimension of 16. 

The important point for our purposes is that the even subalgebra of the Dirac 
algebra Q'^^ is isomorphic to the Pauli algebra Q3 [Hes66|. This isomorphism may 
be constructed by choosing bases 7^ € t/i,3 and o)j G ^3, and defining an invertible 
linear mapping by 

(3-3) oy, e e^3 ^ 7^7t ^ Qt,^ (m e {x,y,z}) . 

These so-called relative spatial vectors 7^7^ satisfy the relations in Eq. ( ^.4| ), since 

KO' (7p7t)' = -7p(7t)'7^ = -(7^)' - 1 
(3.4) (Tf,(T^ ^ (7p7t)(7,.7t) = 7z.(7^7t)7t = -(7i.7t)(7^7t) 

^ - er^a-f, (^, G {x, y, z}, ^ ^ ;/) , 

and hence generate an algebra isomorphic to ^3. As bivectors in Gi,3, however, 
they also generate G13, since 

(3-5) 7^7,. = 7^(7t)^7i. = -(7M7t)(7;.7t) ^ ~(t,,(t^ 

(p, V G {x, y, z}, /i 7^ i^), and similarly 

(3-6) t = 7t7x7y7z = (7x7t)(7y7t)(7z7t) ^ (T>,(Ty(T^ ■ 

Thus 7^7t <-> C/j (m G {^lYi^}) induces an algebra isomorphism as claimed, and 
when the bases are understood we may identify cr^ = 7^7t. 
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The choice of time-like vector 7^ G 9ti,3 in fact determines an inertial frame up 
to spatial rotation, in which the time t and place s of an event e in that frame are 
given by 

(3.7) t + s = e-7( + eA7t = ef^ 

(note that upright case is used for the space-time vector e G 9^1,3) ■ Thus the 
invariant interval between events separated by the space-time vector e is — 
6 7(6= {t + s){t — s) = P — as usual, while the relative velocity between events 
whose space-time velocities are 7^ and v = de/dr is 

ds ds dr f de \ f de v A 7t 



dt dr dt ~ \dT "J \dT 'V v-7t ' 

so that V ■ 7( lies on an affine hyperplane in space-time. 

A great deal of physics can be done in a manifestly Lorentz covariant fashion 
using the Dirac algebra. For example, the electromagnetic field at a given point in 
space-time corresponds to an arbitrary bivector F G called the Faraday 

bivector, and the covariant form of the Lorentz force equation is 

(3.9) mv = qF • V , 

where m is the rest mass, q the charge and v the space-time velocity. (This is 
another example of t he general ru le that, in geometric algebra, the generators of 



motion are bivectors [ DHSvA93 .) The usual frame-depen dent fo rm is recovered 



by splitting the quantities in this equation by 7( as above [ JanSS j; in particular, 
the Faraday bivector splits into an electric and a magnetic field a.s F = E + lB, 
where 

(3.10) E = (F-7j7t and lB = (F A 7^)7^ . 

The space-time reverse will be denoted by a tilde, e.g. in the present case F = — F. 
This is related to the spatial (or Pauli) reverse hy = E — lB = 7^F7^. Both 
operations agree on the Pauli-even subalgebra, but the spatial reverse not Lorentz 
coveriant since it depends on a particular 7^ . 

Returning to our previous discussion of the density operator, we observe that 
the space-time form of the density operator of a single qubit polarized along z can 
be written as 

(3.11) p = i(l-haoi) = i(7t + a7j7t = £.7^ , 

where — 1 < a < 1 is the polarization and 7( determines the local inertial frame. It 
follows that the Lorentz covariant form of the density operator is a time-like vector 
Q G ^1,3- Under a Lorentz boost L = exp(— Acrz/2) G S0(1, 3) along cr^ , therefore, 
the relativistic density operator q transforms to q' = ^i(7t + oi^z)^ = 

(3.12) i(cosh(A)7t-sinh(A)7, + a(cosh(A)7^-sinh(A)7t)) . 

This implies that in the unaccelerated frame (with renormalization by g' ■ 7^), 
^, ^ Q'lt ^ ■ 7t + g' A 7t 
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It follows that the polarization itself transforms as 

, acosh(A) - sinh(A) 
cosn(A) — Q:snin(A) 

If we assume the qubit is at equilibrium with a heat bath, statistical mechanics 
tells us that a = tanh(— /3e/2) where /3 = l/{kBT) is the inverse temperature and 



e S M is the energy difference between the | 0) and 1 1) states | Tol38 |. Then the 
addition formulae for cosh and sinh give 

(3.15) a' = tanh(-/3e/2 - A) , 

so the apparent equilibrium polarization depends on velocity. These results are not 



to be found in the classic treatise on relativistic thermodynamics | Tol34 . 

We will now construct a Lorentz covariant multiparticle theory of qubit systems 
in the simplest possible way, by taking a direct sum of copies of space-time (regarded 
as a vector space, rather than an algebra), one for each of the N qubits, i.e. 

(3.16) e,li (7?,7x^7?,7z^) , 

and considering the associated geometric algebra Gn.sn- Then the even subalgebras 
of different particle spaces p ^ q commute, since (in any given bases) 

(3.17) <<tJ = ilhlj!H = 

for all 11, V ^ {x, y, z}, so that the algebra generated by the even subalgebras is 
isomorphic to a tensor product of these algebras, written as 

(3.18) {gl^r"" « gf"" EE {g^.f' . 

This construction of the tensor product was first used by Clifford as a means of 



studying the tensor products of quaternion algebras [lCli78|; van der Waerden has 
in fact called it a Clifford algebra of the second kind | vdW85 |. As a means of 
justifying the tensor product of nonrelativistic quantum mechanics in terms of the 
underlying geometry of space-time, however, it is a much more recent development 
|DLG93[ 



A key feature of quantum mechanics, which is needed for quantum computers 
to be able to solve problems more efficiently than their classical counterparts, is an 
exponential growth in the dimension of the Hilbert space of a multi-qubit system 
with the number of particles involved. The complex dimension of the Hilbert space 
(i})®^ of an A^-qubit system is in fact 2^, and the space of operators (linear trans- 
formations) on {f))®^ therefore has real dimension 2^^+^. The above construction 
yields a space of "operators" gf^ whose real dimension also grows exponentially, 
but as 2'^^. The extra degrees of freedom are due to the presence of a different unit 
pseudo-scalar l'^ in every particle space. They can easily be removed by multiplying 
through by an idempotent element called the correlator: 
(3.19) C = i(l-/,it2) ... i(i„,i^W) 

This commutes with everything in gf^ and satisfies lPl'^C = C for 1 < p,q < N, 
so that multiplication by it homomorphically maps gf^ onto an ideal gf^/C 
wherein all the unit pseudo-scalars have been identified,n and which therefore has 



^The notation Qf'^/C is justified by the fact tfiat the two-sided principle ideal (C) 
irated by C is isomorphic to the quotient algebra 

gC = 0}. 



generated by C is isomorphic to the quotient algebra 5|^^/kcr(C), where ker(C) = {96 | 
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the correct dimension over M. As a subalgebra, this ideal is in fact isomorphic to 
the algebra of 2^ x 2^ complex matrices, and hence capable of describing all the 
states and transformations of (ensembles of) N qubit systems. In the following, 
we shall generally omit C from our expressions altogether, and use a single unit 
imaginary l as in conventional quantum mechanics. 

On the "even" subalgebra {G^ )'^^ , multiplication by the correlator turns out 
to be an algebra automorphism; this algebra can thus be written as 

(3.20) (^3+)^^ « iOrf'/C « (Gf^'/Cf « 5W(2)«^ , 

where the "+" refers throughout to the subalgebra generated by expressions which 
are invariant under inversion in the origin, and SU{2)^'^ to the algebra generated 
over R by the Kronecker products of matrices in the group SU(2). This subalgebra 
has real dimension 2^^, but is mapped onto a left-ideal of dimension 2^"*"^ by 
right-multiplication with another idempotent which is given by the tensor product 
of those considered earlier, namely 

(3.21) E+ = E\E%---EI , 

where E^. = {lzLa^)/2 for q — 1, . . . , N . Henceforth, the term "even subalgebra" 
will refer to (G^)^'^ (suitably correlated) unless otherwise stated. 

In terms of the usual matrix representation, right-multiplication of an element 
of the even subalgebra ^ by likewise sets all but the first column to zero, so 
that "ifE^ transforms like a "spinor" in under left-multiplication by single 

particle rotations i?' G {G^)'^^ ■ Unlike the single particle case, however, this one 
column does not uniquely determine an element of the even subalgebra (CJ^)®^/C 



What has been proposed instead |DLG93| is to use the fact that E+ "absorbs" 



oi's to distribute copies of the latter across the correlator, converting it to what 
will here be called the directional correlator D, i.e. 

(3.22) ^CE+ = 'irC (^{(T^ f^^cT^ ■■■CT^^^ E+ = ^DE+ , 
where 

(3.23) D = i(l-tV,itV2)i(l-tV,iiV3) ••• i(l-tV,itV^) . 

It can be shown that right-multiplication by D, unlike C, reduces the dimen- 
sionality to 2^+^, thereby permitting the objects in this reduced even subalgebra 
{Q^)"^^ / D to be regarded as spinors, analogous to Gt for '"^ single qubit. In the cor- 
responding left-ideal, K = u^cr^D ^ ■ ■ ■ ^ c^a^ D serves as the unit imaginary, 
since K — —D, but is required to always operate from the right Henceforth, un- 
less otherwise mentioned, we will regard spinors (^3+)^^/D = ((5+)®^/C)/(Cr>) 
as a left-ideal in the C-correlated even subalgebra, drop both C and the super- 
scripts on the l's as above, and use £> as a short-hand for CD = DC. 

In the case of two qubits, for example, the identifications are induced by D are 

100) -iJ^Lcrha^ Lcr} ^ Lcr^ 



(3.24) I 01 ) i^cr^ < > Lcr^Lcr^ —l<t^< >L(T^Lcr^ 

I \ 1 1 ? 1 1 ? 
10) tCTy^ < » terser/ -tcr^ < > i^(Tytcr^ 

ii\ 12^ 12 12^ 12 

II y — LCFy LCTy < — ■ — > LCT^ ^ — — ^ ^^x 
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(where the two columns differ by operation with K) . From this it may be seen that 
any "spinor" in (t/^)®^/£) can be written as 



(3.25) 



(cf. (SLD9§). Alternatively, again using Eq. ( |3.24| ), a unit norm spinor may be 
factorizcd into a product of entities in the correlated and reduced even subalgebra, 
namely <b = R} TPDC, where 

fii = e-^^-zV2 g-.K/2 ^ J. ^ cos(^/2) - sin(^/2) cTyV^K , 
Thus when ^ = tt, the factor T becomes 

(3.27) -<t>2k = {-LCT^^){-LCT^)K = e-(-/2W e-(-/2).-,^ ^ ^ 

so that the arguments of the exponentials involving Lcr^ and Lcr^ in the first two 
factors are shifted by 7r/2 while the total phase is shifted by r = tt. It follows 
that T rotates the first two factors in the planes defined by their conjugate spinors 
[— r2;r*], [— S2;'Si]- Thus on right-multiplying by and expanding in the usual 
basis, we obtain (up to an overall phase) 



(3.28) 



cos(§)e'"/2 
sin(§)e~*"/2 



cos(f)e''^/2 




cos(f)e*^/2 


sin(f)e-''^/2 




sin(f)e-*^/2 



sin(§)e''^/2 
-cos(|)e-^'^/2 



+ 



sin(|)e*^/2 
-cos(|)e-*^/2 



This is known as the Schmidt decomposition |EK95|. It is useful in studying the 
entanglement of bipartite quantum systems, which (in conventional terms) means 
that I ) S cannot be written as a product | tA"^ ) ^ | "0^ ) = I V'"^ ) I ) = \'<P^ "^P"^) 
for any one-particle spinors \ ip^ ) , \ il^'^ ) CI Sj . In fact it is just the singular value 
decomposition in disguise, since (for example) on arranging the entries of a two- 
qubit spinor lip) — [ipi ; ■ • • ; ipi] in a 2 x 2 matrix, we can write 

V'l V's 

1p2 Ipi 



(3.29) 



2Nt 



where V is a 2 x 2 diagonal matrix containing the singular values v^^ > > and 

U, W are unitary matrices with columns u'', w'^, respectively. Since the entries 

of the dyadic products u^^w^f, u^{w^f are exactly the same as the Kronecker 
,1 



matrix products y} ^w}, ^2i^, the equivalence with Eq. (3.28) follows with 



,11 



cos(^/2), V = sin(<^/2), and the Kronecker products of the columns of U 
and W identified with conjugate pairs of single qubit spinors whose relative phases 
are given by exp(±ir/2). 

Clearly a two-qubit spinor is unentangled if and only if v^^ — 1, which is 
equivalent to <j = or T = 1. Thus T describes the entanglement of the qubits, and 
is accordingly called the tangler. The geometric algebra approach clearly provides 
deeper insight into the structure of entanglement than does one based on mechanical 
matrix algebra. In particular, the fact that vEf* is even and reversion-symmetric 
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in the Dirac as well as the Pauli algebra implies that it is the sum of a scalar and 
a four-vector in the two-particle Dirac algebra Q2,e ■ Since Lorentz transformations 
of the spinors cancel, this entity is in fact a Lorentz invariant, and dividing out the 
total phase P as P{^'if )P yields the square of the tangler directly. The availability 
of such powerful methods of manipulating entities in the multiparticle Dirac algebra 
promises to be useful in finding analogs of the Schmidt decomposition for three or 
more qubits. 

4. Quantum Operations on Density Operators 

Quantum computers operate on information stored in the states of quantum 
systems. The systems are usually assumed to be arrays of distinguishable qubits 
(two-state subsystems) , whose basis states | ) and | 1 ) correspond to the binary 
digits and 1, respectively, while the operations are usually taken to be unitary. 
General unitary transformations of the qubits are built up from simpler ones that 
affect only a few qubits at a time, which are called quantum logic gates. The 
representation of these gates in suitable products of Clifford algebras has been 



described in Refs. |SCH98, VlaOl . The goal here will be to show how gates act 
upon spinors in the even subalgebra, and how they can be extended to a wider class 
of nonunitary quantum operations on density operators. 

Given the isomorphism between the algebra of 2^ x 2^ matrices over C and 
Qf^/C relative to a choice of basis in each particle space, it is straightforward 
to interpret matrices in the former as geometric entities in the latter. A matrix 
C/e U(2^), however, does not generally correspond to an entity U in the even sub- 
algebra {g+)'^^/C, so that J7* ^ {g+f'^/D for a general spinor * e {g+)'^'^/D. 
Nevertheless, letting E- = Y[q Et be the idempotent "opposite" to and not- 
ing that this satisfies E+E- — 0, the product of C/* with E+ may be written 
as 

(4.1) U^E+ = {U'^E+ + U^E^)E+ = 2{U^E+)^E+ , 

where the "hat" on [/ denotes its image under inversion in the origin (so that 
= E-), and hence (_)+ is a projection onto the even subalgebra. Because 
{g^)®^ /C and U(2^) are both (2^''^)-dimensional, nothing is lost in this projection! 
Thus we can drop the right-factor of as usual, and define the action of U on 

* e {g+)®^/D as 

(4.2) iJo * = 2((7*£;+)_^ . 

More generally, the usual action of the Pauli matrices on spinors corresponds to 



the following action of the basis vectors on the reduced even subalgebra [DLG93|: 

(4.3) cr^ o * = cr^^oi , t, o * = t^cTz 

The simplest logic gate is the NOT of a single qubit, which operates on the 
computational basis as follows: 

(4.4) iV|0) = |1) ^ -tcTy , N_\l) = |0) ^ 1 

Thus it might appear reasonable to represent the NOT by AT = LCTy E SU(2), but 
when LCTy is applied a superposition (1 — t(Ty)/V2 ^ (| 0) + | 1 ))/\/2, we get (1 -|- 
tcry)/V2 ^ (I ) - 1 1 ) ) /V2 instead of (| ) -f 1 1 ) ) /V2 again. For a single qubit this 
difference is just an overall rotation by tt about oi, but a second qubit can be affected 
by this phase difference between the first qubit's states. Therefore the correct 
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representation of the NOT gate in SU(2) is actually AT = itcTx , which preserves 
this superposition up to an irrelevant overall phase shift: ig;^{ \ 0) + | l))/\/2 

(tcTx) o (1 — tcry)/\/2 = t o (Tx o (1 — t,cry)/\/2 

(4.5) = t o (crx(l - /,cry)oi/\/2) = i. o (-tcTy + l)/V2 

= (-tcry + l)icr,/V2 ^ -z(|0) + |l))/%/2 

More interesting logical operations on the qubits must be able to transform the 
state of one conditional on that of another. The usual way in which this is done is 
via the c-NOT or controlled-NOT gate. As a matrix in SU(4), this is represented 
in the computational basis by 



(4.6) 



2|1 _ 



which makes it clear that this operation NOT's the second qubit whenever the first 
is 1. The corresponding operator in geometric algebra is 



(4.7) 



iV^li = (1 



iCT, 



2 {E\ 



El ta^\ 



This may also be written in exponential form as 



(4.8) 



-/,7r(El(l-<T^^)/2 - 1/4) 



Physical implementations of this operation by e.g. NMR typically expand the ex- 
ponential into a product of relatively simple commuting factors which can be per- 
formed sequentially pSTCOO| . 

Note that since (1 — cr^ )/2 is also an idempotent, AT'^'^ differs from 7V^'^ by 
a swap of the x and z axes for both qubits. This self- inverse operation, called the 
Hadamard transform H, is simply a rotation by tt about the (cTx -|- oi)/-\/2 axis. 
Sandwiching TV^'^ by Hadamards — l{(t^ + (T^)/y/2 to just the second qubit 
gives 



jy2_^2|l_f^2 ^ g-<.7r(Ei/f2(l-oi^)ffV2-l/4) ^ ^-L^iE^El 



1/4) 



(4.9) 

so the c-NOT can also be viewed as a rotated phase shift of the state | 11) by tt. 
The Hadamard gate has the important feature of transforming basis states into 
superpositions thereof; indeed, as an element of the even subalgebra, it actually 
represents the spinor of a uniform superposition directly: 



(4.10) 



H\0] 
H\V, 



(tcr^ 

(tCTz 



2D 
2D 



Thus, by using the relations (4.3), we can show that applying a Hadamard to one 
of two qubits in the state 1 11) followed by a c-NOT gate to the other yields the 
entangled singlet state: N_^^^H} \ 11 ) ^ 



(4.11) ^ 



i ((1 + L)El + (1 - l)E'_ Lrr^) o {{ca,' - Lal)i-La^)) 

i ( ((1 + O [lCT,')) i-LCT^) + ((1 - L) O i-,CT^)) {(tCT^) O (-.^-y^))) 



= i(l — Lcr^)L(T^ — ^{lct^ — L(T^)l(T^ 
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y/^{\10) — |01))/\/2 = ^/^|'^/'_ ). "Quantum" gates like H are not, of 
course, found in conventional boolean logic, and are an essential component of 
all quantum algorithms that are more efhcient than their classical counterparts 
|EJ98| , |CEMM98| . Indeed, the c-NOT gate together with general single qubit 



rotations are known to generate SU(2^), and hence are universal for quantum 
logic [1BBC+95[| . 



It turns out that unitary transformations are not the most general sort of 
operation that can be applied to a quantum system. Most such quantum operations, 
however, produce a statistical outcome, and the ensemble of possible outcomes 



must be described by a density operator. The previous definition (Eq. (2.20)) of 
the density operator of an ensemble of identical and noninteracting qubits may be 
extended to an ensemble of multi-qubit systems as follows: 



(4.12) p = <-> P = (*D)-E+(*D)~ = ^E+^C 

Suppressing the correlator C as usual, p may also be expressed in diagonal form as 

(4.13) P - R(j:ll,pkE^k))R^ = ElloPkrkri. 

where R G Qf^/C corresponds to a unitary matrix R G U(2^) (in the usual cr^ 
coordinate system), and < pfc < 1 are the eigenvalues of p. The idempotents 
E^,) are given by Ug ^ \ xl ■ ■ ■ ){xl ■ ■ ■ X^f I where el = 1 - 2x1 with xl 
equal to the q-th bit in the binary expansion of fc G {0, . . . , 2^ — 1}. It follows that 
\ Pk) ^ fk = RE(^k) for pk > are the spinors of the (unique, if pk ^ pe ^ I) 
minimal ensemble that realizes p, which therefore describes a pure state if and only 
if it has rank 1 as an operator. 



Note that by Eq. (4.2), the density operator transforms under unitary opera- 
tions as 



(4.14) p ^ ((7o *)£;+((7o *)~ = U^E+^U^ = UpU'' 
Similarly, the ensemble-average expectation value of any observable O — G 



gf'/c is 



(4.15) (^|0|^) <-> 2" {E+<IfO^E+), = 2''{O^E+^), 2'' (Op), 



just as shown in Eq. (2.17) for single qubit ensembles. In contrast to the case of 
a single qubit, however, the geometric interpretation of these observables is not 
straightforward. While one can certainly express p as a finite ensemble average 
X)fcPfc*fc-^+*fc (where the pk > are probabilities with ^i^Pk = 1), this decom- 
position is highly nonunique. The minimal ensemble obtained by diagonalization, 
on the other hand, will generally include entangled spinors , for which the ex- 
pectation value (Orj,rJ.)p cannot be expressed as a product of inner products of 
the factors of O = ■ ■ ■ with the polarization vectors of the individual qubits 
(indeed, O itself need not be factorizable!). 

The best one can do is to expand p in the product operator basis consisting of 
all 2^^ products of the basis vectors cr^, i.e. 

(4-16) p= ^ Pki^-t.^ , 

Aii,...,A'"6{0,x,y,z} 

where y0^i...^« S M and ctq = 1 for notational convenience. The utility of this basis is 
most simply demonstrated via a concrete example, namely NMR spectroscopy. Here 
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one is given a liquid sample consisting of identical molecules whose nuclear spins 
are chemically distinguishable, and hence constitutes an ensemble of multi-qubit 
systems (see |CPH98| , [HCSTO0| , [HSTCO0| and references therein). The energy 



of interaction between the spins and an external magnetic field along z is given by an 
observable called the Zeeman Hamiltonian, Z = [lo^ct^ + • • • + uj^ (t^)/2, where u'^ 
is the energy difference between the | ) and 1 1 ) states of the q-th spin in the field. 
In thermal equilibrium at room temperatures, the polarization of the spins relative 
to the strongest available fields is typically a ^ 10^^, and the density operator of 
the ensemble is essentially p^^ = {l + a{(T^ + ■ ■ ■ + (t^)). Via a suitable pulse of 
radio-frequency radiation, this may be rotated to p^q = 2^^(1 + a{(T^ + • • • + cr^)), 
which evolves under the interaction with the field as 

(4.17) e^^^Veqe'^* = 2'^ {l + a {cos {ujH)(t^ - sin{ujH)(T^^ +■■ ■ 

■■■+ cos(c^^t)of - sin(L^^t)of )) . 

Thus on measuring the total magnetization Mx along the x axis, O = ^{cr^ + 
• • • + (T^) (where 7 is the nuclear gyromagnetic ratio), one obtains the sum of the 
projections of the rotating magnetization vectors of the spins along the x-axis, i.e. 

(4.18) M^(t) = Q!7(cos(cj4) H h cos(tj^i)) , 

whose Fourier transform reveals the contribution from each spin. The way in which 
the factors of product operators transform like vectors under rotations accounts in 
large part for the computational utility of the product operator basis. Of course, 
unless it is a natural part of the problem at hand (as in NMR), one is better off 
not chosing a basis at all! 

A normal quantum operation is a linear transformation of the density operator 



that may be written in operator sum form as [Kra83] 

(4.19) p ^ n{p) EE EkQkpQl , 

where the Kraus operators Qf. G Gf^/C satisfy J2k QkQk = 1- The term "normal" 
here^ refers to the fact that such an operation preserves the scalar part of p, since 

(4.20) {n{p))^ - Y.k{QkPQl)^ = (pEfcQlQfe)^ = 2"^ • 

It is also easily seen that such quantum operations are positive, in that they preserve 
the positive-definiteness of p; in fact, these operations have a yet stronger property 
known as complete positivity, meaning that if the qubits to which f2 applies are 
embedded in a larger system, then applying Q to just those qubits preserves the 
positive-definiteness of the larger system's density operator. That this is a nontrivial 
extension of positivity is shown by the two-qubit "partial transpose" operator flj, 
which carries cr^ i— s- —cTy but leaves all the other operator factors unchanged; this 
is clearly positive on density operators not involving the second qubit, but acts on 



the density operator of the singlet state (Eq. (4.11)) as 



EE lia^ + - al - a^)E\El{cT^ + - crl - cr^) 
(4.21) =\{l~a^a^-a^a^-CT^crl) 

which has eigenvalues [1/2, 1/2, 1/2, —1/2]. 



' We prefer to avoid the more common but clumsy and matrix-bound term "trace-preserving" . 
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A quantum operation Q is called unital if it preserves the identity itself, i.e. 
f2(l) = 1, or equivalently, X^/cQfcQfc ^ ^- Perhaps the most important example 
of a normal unital operation is found in the contractiom by a single qubit q E 



{1, . . . , N}, which may be written in operator sum form as |SCH98| 



(4.22) 2 (p)« = ElpEl + ElpEl + a^{E\pE\ + ElpEDtr^ 

This may also be expressed by dropping all terms in the product operator expansion 
of p depending on g, and multiplying the remaining terms by a factor of 2. Note 
that, while (_)^ is normal and unital, this factor means that the contraction itself 
is neither. The factor is nevertheless required if the result is to be interpreted as 
a density operator for the remaining qubits, since the contraction by the second 
qubit of the above singlet state is )^ — 1/4 (not 1/2). 

This example also illustrates an important way in which general quantum 
operations are realized in practice, despite the fact that the universe as a whole 
evolves unitarily. As shown previously, the superposition state with spinor = 
(1 — tcTy )/-\/2 is converted into the singlet state with density operator ■0_ by let- 
ting it interact with a second qubit so as to effect the c-NOT operation 7V^'^. The 
contraction then corresponds to "discarding" the second qubit (i.e. ensuring that it 
does not further interact with the first and hence can be ignored) , which yields the 
density operator 1/2 of the totally mixed state for the first qubit. Since the basis 
states E]. are unaffected by 7V^'^, the net quantum operation on the first qubit 
corresponds to what is known in quantum communications theory as the phase 
damping channel 

(4.23) p^ {1 - p)p + pE+pE+ + pE^pE^ 

with damping parameter p — 1. Phase damping is also known as T2 relaxation in 
NMR, and as decoherence in quantum information processing; it is widely believed 
to be the dominant mechanism by which classical statistical mechanics arises from 



the underlying unitary dynamics |GJK+96 



To illustrate the utility of geometric algebra in the study of general quan- 
tum operations, an eigenvalue characterization of normal, unital, one-bit quantum 
operations f2 will now be derived. This characterization was originally given by 
Fujiwara & Algoet | FA99| , although the derivation he re paral lels that more re- 



cently obtained using matrix methods King & Ruskai |KROO| . This derivation 
will regard the Kraus operators Q). £ G3 as "complex quaternions" Ak + iBk with 
Ak , Bk G , and consider the action of an arbitrary operation O on the scalar 
and vector parts of p = (1 + r)/2 separately. 
First, the action on 1 is 

^^(1) =T.ki^k + '^Bu){Au + iBt,)^ 
(4.24) , . . , , 

= Efe AkAu + BkBk )+LY.k\ BkAk - AkBk) 



The first summation is symmetric with respect to spatial reversion and inversion, 
i.e. scalar, while the second (excluding the t) is reversion antisymmetric but in- 
version symmetric, i.e. a bivector. Writing A = a + La and B = P + tb, so that 



^Otherwise known as the "partial trace" 
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AkAk + BkBk = a\ + lla/clP + 01 + ||b/c|P, we may further expand the second 
summation as follows: 

Sfc [{Pk + tbk){ak - Lttk) - (ttfc + Lak)(f3k - tbk)] 

= J2k i'^kPk + i-iakbk - /JfeOfe) + bfeOfc) - 

(4.25) L 
{akPk + i^ihak ~ akbk) + akbk) 

= Y,k [ '^'^{c^kbk - PkCik) - 2ak A bk ] 
Thus is unital if and only if 

(4.26) Y.k i.'^kbk " PkO-k) = Y.k "■k X bfc , 
(where "x" denotes the cross product) and 

(4.27) Y.kH+wr+Pi + \\bkr) = 1 . 

Similarly, the action on r e 9^3 is 

^{r) ^Y.ki^k + iBk)r{Ak-LBk) 

(4.28) / . . N / , 
= Y.k ^kvAk + BkrBk) +Lj:ki BkvAk - AkxBt 



The first summation is over different dilation/rotations of r; the second summation 
(excluding the t) is reversion and inversion antisymmetric, i.e. a trivector, and may 
be further expanded as above: 

Y,k [ iPk + tbk) r {uk - iak) - {ak + tak) r {(3k - ibk) ] 

= Hk {akPkT + L{akbkr - Pkrak) +bkrak) - 

(4.29) ^ 
{akPkr + tiPkO-kr - akvbk) + akrbk 

= J2k [ '^i-'^kbk ■ r - 2Lj3kak • r + 26^ A r A ] 
Multiplying through by — 1./2 converts this to 

(4.30) J2k [o^kbk ■ r - Pkak ■ r - {L{ak A bk)) ■ r] , 
which vanishes if and only if 

(4.31) Y^ki^kbk - PkCtk) ■ r = Y.k(^k ^ dk) ■ r . 
If Q is normal, this must be true for all r, which is equivalent to 

(4.32) J2k ( ^kbk - PkCLk ) = J2k^k^ ak ■ 



A comparison with Eq. ( [f.2q ) shows further that is both unital and normal if 
and only if J^k ^kbk = J2k Pk^k and Y.k ^k x bk = 0. 

If we regard a normal f2 as an affine transformation, i.e. 

(4.33) ^ikil + r)) = i(l + < + E^K-^)^A.) 

{Sfj_,t e 9^3), we see from the derivation leading up to ( 4.2(^ ) that 
t = 2j2ki(^kbk - PkO-k - ak ^ bk) 

(4 34) 

^ ^ =Aj:ki^kbk-Pkak) = ^Ekbk^ak (by (4.32)) . 

Similarly, the above vectors = {fl{cr^))^ are 

(4.35) Sf, = J2k(^k(Tf,Ak + Bk(T^Bk) ine{x,y,z}), 
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i.e. a sum of independent dilation/rotations of each basis vector. 

It follows that normal and unital quantum operations may be characterized 
by finding conditions for the linear map r i— s- {fl(r))-^ — ^{r) to be written as a 
sum of dilation/rotations of r. To this end, we expand AkvAk as 

(4.36) ^ ^ ^ ^ ^ ^ 

= Q^r + 2akr x ak + 2(r • ak)ak - ||afc|| r , 



with a similar expansion for B^rBk ■ Thus on assuming that U, is diagonal, i.e. 
Apcr^ = fi(cr^), we get 

(4.37) 



'^M EfeK + - llafell^ - llbfell^) + 2a-^ X EfcK«fe + Pkbk) 
+ 2l]fc ((cr^ • ak)ak + (oy, • bfe)bfe) . 
Dotting both sides by cr^ now yields 

A,, = Efc K + - lla.lP - llbfelP + 2 {cT^ ■ a^f + 2 {a,, ■ bkf) 
(4.38) = 1 - 2Efc dlafcll' + llbfcll' - (t^ • akf - {cr^ ■ bkf) 

= 1 + 2 Efc (('^A. A afe)2 + [a^ A 6,-)') , 
so we have simple expressions for the eigenvalues. Now consider the vector obtained 
from the first line of this equation, i.e. E;^G{xyz} -^^0)1 



(4.39) 

where 
(4.40) 



E^ Ea 
Po Efc (a? 

Po = 



Pi 



-Pi) 



Qfcll^ - llbfell^ + 2 (cr^ • Ofe)^ + 2 (cr^ • bfc)^ 
E^Pm Efc (K< • afe)^ + {<T^, ■ bkf) , 



Px 
Pz 



Since the coefficients of the p's are nonnegative and sum to 1 by ( 4.27 ), this shows 
that the vector A^<t^ lies within the tetrahedron (pQ,p^, Py,p^) , which is the 
condition on the eigenvalues found by Fujiwara and Algoet [ FA99 as well as by 



King and Ruskai [KROO 



It is also known that an arbitrary linear map VL has an operator sum repre- 
sentation if and only if it is completely positive |Sch96], so the above can also 
be viewed as a characterization of complete positivity for normal and unital maps 
of a single qubit's density operator. Finally, it is worth stressing once again that, 
because of the isomorphisms which exist between the Pauli algebra and the even 
subalgebra of the Dirac algebra, every step of the above derivation carries with it 
a natural interpretation in space-time, and is in fact even easier to carry out when 
the full power of the Dirac algebra is used. 

In conclusion, it is hoped that the forgoing has given the reader a taste of 
the new insights which geometric algebra can provide into quantum information 
processing — and an appetite for more! 
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